On ambiguities of sign determination of the S-matrix from energy levels in a finite 

box. 



(N 

o 

(N 

> 
O 



43- 



>: 

in, 
oo 
ov 



(N 



E. Oset 

Departamento de Fisica Teorica and IFIC, Centra Mixto Universidad de Valencia- C SIC, 
Institutos de Investigation de Paterna, Aptdo. 22085, 46071 Valencia, Spain. 
(Dated: November 19, 2012) 

In a recent paper the authors make a study on the determination of the S-matrix elements for 
scattering of particles in the infinite volume from the energy levels in the finite box for the case of 
multiple channels. The study is done with a toy model in 1+1 dimension and the authors find that 
there is some ambiguity in the sign of nondiagonal matrix elements, casting doubts on whether the 
needed observables in the infinite volume can be obtained from the energy levels of the box. In 
this paper I present an easy derivation, confirming the ambiguity of the sign and argue that this, 
however, does not put restrictions in the determination of observables. 

PACS numbers: 



I. INTRODUCTION 

The determination of scattering amplitudes and 
hadron spectra is one of the challenging tasks of Lat- 
tice QCD and many efforts are being devoted to this 
problem [l- 26J . For the case of one channel scattering, 
and resonances decayin g in jus t one channels, Liischer's 
approach is often used [23, UM ■ The method allows to 
reproduce the phase shifts from the discrete energy lev- 
els in the box. This method has been recently simplified 
and improved in Ref. [29( by keeping the full relativis- 
ts two body propagator (Liischer's approach keeps the 
imaginary part of this propagator exactly but makes ap- 
proximations on the real part) and extending the method 
to two or more coupled channels. The method has also 
been applied in Ref. (30j to obtain finite volume results 
from the Jiilich model for meson baryon interaction and 
in Ref. [3l| to determine the strategy to find the two poles 
of the A(1405) from lattice QCD simulations. Further ap- 
plications and references to works done along these lines 
can be seen in 32]. 

The extension of Liischer's approach to coupled chan- 
nels has been addressed in several works [29|, [H, [H, [H, 
[55j . In the work of Ref. [29[ , the inverse problem of get- 
ting phase shifts and resonances from lattice results using 
two channels was addressed, paying special attention to 
the evaluation of errors and the precision needed on the 
lattice results to obtain phase shifts and resonance prop- 
erties with a desired accuracy. 

In a recent paper (3?| a coupled channel study of this 
inverse problem is done with two channels, using a toy 
model in 1+1 dimension and it is concluded that the de- 
termination of the S-matrix in the case of T-invariance 
(we only consider this case) has an ambiguity in the sign 
of the nondiagonal matrix element. Due to this, doubts 
are cast that in a realistic case this does not pose prob- 
lems in the determination of needed observables. 

In the present paper we present a derivation of this 
inverse problem in the realistic case and conclude that 
indeed this ambiguity arises. However, we also argue 
that this ambiguity, tied to an arbitrary phase in the 



wave function of one channel with respect to the other, 
has no repercussion in the determination of observables. 



II. FORMALISM 

In the chiral unitary approach the scattering matrix in 
coupled channels is given by the Bethe-Salpeter equation 
in its factorized form [38|. We assume just s- waves for 
simplicity (see generalization to other waves in [3{|) and 
we have 



T=[l- VG]~ l V = [V~ 



G\ 



(1) 



where V is the matrix for the transition potentials be- 
tween the channels and G is a diagonal matrix with the 
i th element, Gj, given by the loop function of two prop- 
agators (we shall use two mesons), which is defined as 



Gi = i 



1 



d 4 p 

(2tt) 4 {P-p) 2 -Mf+ie 



1 



, (2) 



where rra and Mi are the masses of the pair of mesons 
and P the four-momentum of the global meson-baryon 
system. 

The loop function in Eq. © needs to be regularized 
and this can be accomplished either with dimensional 
regularization or with a three-momentum cutoff. The 
equivalence of both methods was shown in Refs. (40l.l4l|. 

In the cut off method a cutoff in three momentum is 
used once the p° integration is analytically performed 
[IH, and one gets 



G; = 



d 3 p 



1 



+ oj 2 (p) 



|p1<Pn 



(2tt) 3 2 Wl (p)u; 2 (p) £ 2 - + tu 2 {p)Y + ie ' 



uiaip) = \J m l,2 +P 2 > ( 3 ) 
with mi , to 2 corresponding to TOj and Mi of Eq. @ . 
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When one wants to obtain the energy levels in the fi- 
nite box, instead of integrating over the energy states of 
the continuum, with p being a continuous variable as in 
Eq. ([3]), one must sum over the discrete momenta allowed 
in a finite box of side L with periodic boundary condi- 
tions. We then have to replace G by G = diag (Gx, G2) 
(in two channels), where 



-, |p|<Pm £ 

p 



1 



0Jl(p) + L0 2 {P) 



2wi (p)w 2 (p) # 2 - (Wi(p) + L0 2 {p)f 



p 



L 



n S 



(4) 



This is the procedure followed in Ref. [29(. The 
eigenenergies of the box correspond to energies that pro- 
duce poles in the T matrix, Eq. ([T]), which in the finite 
box correspond to zeros of the determinant of 1 — VG, 



det(l - VG) = 0. 



(5) 



For the case of two coupled channels Eq. ((SJ) can be writ- 
ten as 



det(l - VG) 



: 1 - VnG, 

+ (VnV 2 2 - 
0. 



- V 22 G 2 
^12)^1^2 



(6) 



One can already see there that V12 appears squared, 
hence, a change of sign in it will not change the spectra 
of levels in the box. We shall see that this is also the case 
in the inverse problem and we can only determine T^ 2 - 
However, this will not prevent us from determining the 
three scattering magnitudes in this case, the two phase 
shifts and the inelasticity. 



III. THE INVERSE PROBLEM 

The inverse problem of obtaining Tij from the energy 
levels in the box is most efficiently written in terms of 
SG defined as 



5G = G — G, 



(7) 



the magnitude used in [42j], which is finite in the limit of 
Pmax 00, the limit taken in [3ll l42l|. and proportional 
to the Liischer function [2^, l43j |. 

Let us start from Eq. |T|) that gives the T matrix in the 
infinite volume and write the correspondent scattering 
matrix in the finite volume, T: 



T(E) = [G(E)-G(E)Y 



(9) 



and this is the formulation of Liischer's formula in [29| . 

Coming back to the multichannel problem and using 
Eqs. ([T]) and ©, as done in [3l|, we get 



Hence, 



f- 1 = T- 1 - SG = T -1 [l -T8G], (10) 



T =[1 — T8G]~ 1 T. (11) 



which allows us to get T directly in term of T, without 
going through an intermediate potential. One can note 
that this formula is like the one of Eq. ([T]), or Eq. ([5]) 
for f, substituting V — > T and G —>■ SG. Hence, the 
condition to obtain the energy levels in the box, det(T) = 
0, leads to the analogous secular equation of Eq. ^ in 
terms of T and SG substituting V and G, respectively, 

(1 - T u SG n )(l - T 22 SG 2 ) - T? 2 5G 1 6G 2 = 0, (12) 

or equivalently, 



Tn SGu + T 22 SG 2 - {T 1X T 22 - Tl 2 )SG x SG 2 = 1. (13) 

It is clear that with just one energy eigenvalue for a 
given L, Eq. (|13p cannot provide the full matrix. 
In [29( different methods were used to get T from the 
eigenenergies of the box. The simplest conceptually is to 
take a certain energy from three levels, which correspond 
to three different values of L, and then determine the 
three values of T^. Actually, what one determines is 
Tn, T22 and T^ 2 - So, we can see in a realistic case that 
we only obtain T^ 2 - However, this is not a problem to 
determine the observables. This indetermination should 
be there because a change of relative sign of the wave 
function for the states 1 and 2 leaves Tn, T 22 unchanged 
but it changes the sign of T12, and the physics cannot 
depend on this sign. To make this more explicit we write 
explicitly the S-matrix in terms of the observables in the 
next subsection. 



A. Phase shifts and Inelasticities: 

In order to obtain the phase shifts and inelasticities we 
use the two-channel S matrix yjl 



f = [V- 1 - G]~ 



(8) 



Note that for the case of one channel the poles of G 
provide V~ x — G = and, thus, V^ 1 = G for the eigenen- 
ergies of the box. Then we can recast Eq. ([T]) as 



S = 



qjiti j(l _ ^2)1/2 e i(S 1+ S 2 ) 



(14) 



where S±, S 2 are the phase shifts for the 1 and 2 channels 
and r\ is the inelasticity. The elements in the S matrix 
are related to our amplitudes via: 
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Til 

T22 



= -^?(5n-l) 



2ip 2 ' 



T 2 1 



22 — 1 

87T-y/s 



(15) 



It is interesting to recall that the first two equations 
allow us to determine 8\ , 82 and i] while the third equa- 
tion allows us to determine 8\ + 82 and rj. The S-matrix 
should be unitary and one can see that this is the case 
even if we change the sign of Si 2. In fact the choice of 
the positive sign in the square root of (1 — rj 2 ) is a matter 
of convention. One may argue that if we change the sign 
of Si 2 and use Eq. ijlgj). one would obtain e l{Sl+s ^ with 
opposite sign to what one would get from using the same 
matrix element before changing the sign. This reflects 
the arbitrariness of 7r in the phase shifts. In view of this 
arbitrariness in the sign, the quantities that should be 
used are | X\ 2 1 2 , which determines (1 — rj 2 ) and then T 2 2 
which determines (1 — rj 2 )e 2l ( Sl+6 '* > and allows one to get 
81 +82, independently of the sign of T\2- 



IV. CONCLUSIONS 

In this paper we have done a derivation of the inverse 
method to get the scattering matrix from the energy lev- 
els of the system in a finite box. We observe, in agreement 
with the findings with the toy model of (37| , that for the 
case of two channels studied in [si! , the sign of the off 
diagonal matrix element T\2 is not defined. Yet, we could 
see that this had not repercussion in the determination 
of the observables 5\, 82 and rj, although given the ambi- 
guity on this sign, tied to an arbitrary relative phase in 
the wave functions of the states, a particular algorithm 
must be taken to determine the phase shifts from Ti 2 . 

The two channel system that we had in mind was the 
one of two physical states that couple, say tttt and KK, 
studied for instance in [29( . One may think of other cou- 
pled channels systems, like one physical state with two 
coupled partial waves, where signs and interference are 



important in angular distributions (the authors of [33] 
might have such and idea going beyond the toy model 
used). An example of this is the deuteron or dineutron 
system in the presence of a tensor force [45]. The issue 
of partial wave mixing is an interesting one in finite vol- 
ume because a square box breaks rotational invariance 
and this leads invariably to par tial wave mixing. So, the 
problem is well documented [28[ . Closer to the case of ex- 
plicit L mixing caused by particular external forces is the 
case of the L mixing in the moving frame. Here one still 
has a central potential but the imposition of the bound- 
ary conditions in a frame where the total momentum of 
the system is not zero leads to partial wave mixing. This 
problem has also received much attention [3^, I4614571 ]. 
The exposition and solution of the problem is made in 
a relatively simple and pedagogical way in [391 ] , where L 
mixing and physical coupled channels are considered si- 
multaneously. There one can see that the equations and 
strategies that allow one to obtain the different partial 
ways and their relative signs are far more involved than 
the simple equations used in the present work or in [33]. 
The result of these works is that the inverse problem has 
solution without a mbig uities and realistic cases are even 
solved explicitly in [391 ]. 

Whether in some finite volume topic ambiguities ap- 
pear in the inverse problem is an open issue. What is 
clear is that the ambiguity found in the present work 
dealing with just one partial wave, or in [371 ] in the 1+1 
toy model, where partial waves cannot be defined, poses 
no problem for any of the works done with finite volume 
so far. 
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